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Abstract 

Using the standard saddle-point method, we find an exphcit relation for the large- A'^ 
limit of the free energy of an arbitrary generalized 2D Yang-Mills theory in the weak 
[A < Ac) region. In the strong {A > Ac) region, we investigate carefully the specific 
fourth Casimir theory, and show that the ordinary integral equation of the density 
function is not adequate to find the solution. There exist, however, another equation 
which restricts the parameters. So one can find the free energy in strong region and 
show that the theory has a third order phase transition. 
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1 Introduction 



The pure 2D Yang-Mills (tr(F^)) theory (YM2) defined on compact Riemann surfaces, is 
characterized by its invariance under area preserving diffeomorphism and the fact that there 
are no propagating degrees of freedom. This theory is not unique in this sense, and it is 
possible to generalize it without losing these properties. 

In an equivalent formulation of YM2 theory, one can use iii{BF) + tr(i?^) as the La- 
grangian, where B is an auxiliary pseudo-scalar field in the adjoint representation of the 
gauge group. Path integration over the field B leaves an effective Lagrangian of the form 
tr(F^) [1]. Now, a generafized 2D Yang-Mills theory (gYM2) is a theory with above La- 
grangian, in which the term ti^B^) is replaced by an arbitrary class function A(S) . 

In [2], where these theories were first introduced, the partition function of gYM2 have 
been obtained by considering its action as a perturbation of the topological theory at zero 
area. In [3], these theories have been coupled to fermions, to obtain the generalized QCD2. 
Generafizing the Migdal's suggestion about the local factor of plaquettes [4], the authors 
of [5] have found the partition function and the expectation values of the Wilson loops for 
gYM2's. In [6], these quantities and also the generating functional of the field strength have 
been calculated for arbitrary two dimensional orientable and nonorientable surfaces, using 
the standard path integral method. 

One of the important features of YM2, and also gYM2's, is their behaviour in the case 
of large gauge groups, e.g., the large N behaviour of SU{N) (or U{N)) gauge theories. On 
one hand, it is interesting to study the relation between these theories at large and the 
string theory. Such investigations began in [7] and [8] for YM2 and in [5] for gYM2. In 
these papers, it was shown that the coefficients of 1/N expansion of the partition function of 
SU (N) gauge theories are determined by a sum over maps from a two-dimensional surface 
onto the two-dimensional target space. 

Another thing of interest, also related to the above mentioned one, is to study the large 
A^ behaviour of the free energy of these theories. This is done by replacing the sum over 
the irreducible representations of SU{N) (or U{N)), appearing in the expressions of parti- 
tion function and Green functions, by a path integral over continuous Young tableaux, and 
calculating the area- dependence of the physical quantities from the saddle-point configura- 
tion. In [9], It was shown that the free energy of U{N) YM2 on a sphere with surface area 
A < Ac = n"^ has a logarithmic behaviour. Later, the authors of [10] calculated the free 
energy for surface area A > tt^, and showed that the YM2 has a third order phase transition 
at the critical surface area. Such kind of calculation for gYM2 began in [11], and the authors 
claimed that these theories have a rich phase structure. They also discussed the general 
behaviour of the density function pior A< A^ and brought some comments for the A > Ac 
region. 
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In this paper, we first complete the calculation of the free energy of U{N) gYM2 for 
A < Ac region for an arbitrary gYM2 on sphere and derive an exact solution for the density 
function. But the main part of the paper is about the A > Ac region of a specific fourth 
Casimir model. It is seen that here, unlike the YM2 model, the density function p has two 
maxima, for A < Ac- So some new features arise in the phase transition, which must be 
handled more carefully. At the end we show, however, that for this particular model there 
exist a third order phase transition, just as in the case of YM2. This is another similarity 
between these two theories. 



2 Large- A/^ behaviour of gYM2 at A < Ac 

The partition function of the gYM2 on a sphere is [5,6] 

Z^j:4e-^^^^\ (1) 

r 

where r's label the irreducible representations of the gauge group, dr is the dimension of the 
r'th representation, A is the area of the sphere and A(r) is : 

p 



Mr) = El^.C,ir), (2) 

k=l 



in which Ck is the /c'th Casimir of group, and a^'s are arbitrary constants. Now consider the 
gauge group U{N) and parametrize its representation by rii > 712 > • • • > un, where rii is 
the length of the i'th row of the Young tableau. It is found that [12] 

n (1+^) 

l<i<j<N J ' 

N 

Cfc = E[K + ^-0'-(^-Ol- (3) 

i=l 

To make the partition function (1) convergent, it is necessary that p in eq.(2) be even and 
Gp > 0. 

Now, following [9], we can write the partition function (1) at large A^, as a path integral 
over continuous parameters. We introduce the continuous function : 

(f)[x) = -n{x) - 1 + X, (4) 

where 

<x ■.= i/N <1 and n{x) := rii/N. (5) 
The partition function (1) then becomes 

f Y[ #(x)e^['^(^)], (6) 

0<a;<l 
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where 

S{<f>) = N'{-A r dxG[<t>{x)] + r dx r dy \og\<P{x) - <f>{y)\ }, (7) 
Jo Jo Jo 

apart from an unimportant constant, and : 

= E(-l)'a.0^ (8) 

k=l 

As — > oo, the partition function (7) is determined by the configuration which maximizes 
S. The saddle point equation for S is : 

where: 

g{<P) = ^G'{<P), (10) 
and P indicates the principal value of the integral. Introducing the density 

the eq.(9) is reduced to: 

.W-Pf^^. (12) 

Jb Z — A 

along with the normalization condition 

r p{X)d\=\. (13) 

J b 

The condition ni > n2 > • • • > un imposes the following condition on the density p(A) : 

p(A) < 1. (14) 
To solve eq.(12), we define the function H{z) in complex 2;-plane [13] , 

ifw-r^. (16) 

Jb Z — A 

This function is analytic on the complex plane exept for a cut at [b, a] . There , one has 

H{z±ie) ^ g{z)^i7rp{z) b < z < a. (16) 

H is found to be [10,14] : 

H(z) = ±J(z-a)iz-,,)i , (IT) 

27^^' Jc (^z - X)^{X- a){X-b) 
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where c is a contour encircling the cut [b, a] , and excluding z and g{X) is defined through 
eq.(lO). Deforming c to a contour around the point z and the contour (a contour at the 
infinity ), one finds: 

H{z)=g{z)-J(z-a)(z-b) V (2n - l)!!(2g - 1)!! n^g^m (n+,+m+l)(o^ ^^g^ 

where the n-th derivative of g. It follows from (15) and (13) that H{z) behaves like 

1/z at z ^ oo, or (z — a)^^/'^{z — h)^^^"^ H{z) behaves like l/z"^ at z — > oo. Using this, one 
can expand {z — a)~^/^(z — b)~^^'^H{z), equate the coefficients of 1/z and 1/z^ equal to 
and 1, respectively, and arrive at 

£ P!i^im^„«,,,(n.,)(0) = 0, (19) 



and 



E (^-^f^^^-^^;,a"6V-^--^)(0)^l. (20) 



These equations determine a and b. Prom (16), p is determined : 



V(« -.)(.- 6) ~ (2n-l)!.(2,-l)N 

TT „,t?.„2"+'n!<,!(m + n + , + l)! ^ 

Defining the free energy as 

F - -^In^, (22) 

one has ^ 

F'(A) = / ^[^(x)] = G(X)p(X). (23) 

Expanding H{z) for ^ — > oo, one can find the integrals 

TdA AV(A), (24) 

by using the eqs.(15) and (18). Knowing these, one can calculate F'{A) through (23). 

It is also worth mentioning that if G is an even function, then b — —a, p is an even 
function, and equations (17), (18), (20) and (21) become 

H{z) ^ -j-V^W / , , (17)' 

^ ' 2m Jc (z- A)VA2-a2' ^ ' 

H(z) = giz) - ^/^^ V ; ^" ,, a^''^''ff('"+^+^)(0), (18)' 

^ ' ' „^o2"n!(2n + g + l)! ^ ^ ^' ^ ^ 
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and 



— z 



2 oo 



(2n- 1)!! 



(21)' 



TT „^o2"n!(2n + g + l)! 
As the simplest example of gYM2, consider = 0^. Using the above relations, one finds 



94 ^2 



TT ' 2 



— z^ 



4 



(25) 



F\A) 



1 

4A' 



The density p has a minimum at 2; = and two maxima at z^^ — Now as 

p(4^^) = V2a3A/7r, if A > Ac = 277r''/256 , then the condition p < 1 is violated. So the 
solution (25) is valid only in the region A < A^.. 

It is worth noting that the function F'{A) for G{(j)) = (j)^ can be simply found by rcscaling 
the field by A^/^ in eq.(7) (as noted in [10] ) , after which the only ^-dependent term of 
S is logAV^ and therefore F'{A) = ^. 



3 The G^cj)) = 0^ model in the region A > Ac 

In this region, the solution (25) is not correct any more. As p{z) for A < A^ has two 
symmetric maxima, we use the following ansatz for p : 



1, ze [-b, -c] U[c, b] =: L' 



[ ps{z), z e [-a, -b] U[-c, c] \J[b, a] -. L. 
Putting this in (7), we have 

/"a ra 
dz Ps{z)z^ + dz d'wps{z)ps{'w)log\z — 
-a J —a J —a 



(26) 



w\ 



(27) 



To maximize this along with the condition (13), we introduce the Lagrange multiplier p, and 
the functional 5' : 



S^S + N^p 



f 

J— a 



dz ps{z) - 1 



(28) 



The variation of S is 



5S ^ N"^ [ dz -Az'^ + 2f dw ps{w)\og\z - w\ + p 

J —a L J —a 



Sps+N' 



ra 

/ dz ps{z) - 1 

J —a 



5p, 
(29) 



where 



5ps{z) 



0, 



zeL' 



Sps{z), zeL. 



(30) 
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Equating 5S to zero, we have 

A 



— 2;^-y dw ps{w)\og\z - w\ ^ z e L, 



and 



ra 

/ dz ps{z) = 1. 

J— a 

Differentiating (31) with respect to z, we have 



(31) 
(32) 

(33) 

z — w 

This equation is, however, not equivalent to (31), since L is a disconnected set. In fact, (31) 
is equivalent to (33) and the following equation 

2 



2Az^ 



J —a 



dw 



Ps[W) 



z e L. 



J -a C — W 



(34) 



which is difference of eq.(31) a.t z = b and z = c. The above equation can also be written as 



r 



dz 



2Az' 



■f 

J —a 



dw 



z — w 



= 0. 



(35) 



So our task is to solve the eq.(33), along with the conditions (32) and (35). To do so, we use 
the procedure of the previous section: 



HJz) 



pTdw ^ 

J -a Z — W 



dw h log h log - 

L z — w z + c z — b 



:= Hs{z) + log^— + log 

z + c z — b 



(36) 



where Hg has a three cut singularity a,t z E L. Similar to eq. (17), the solution of Hg is [14] 



, (37) 



where 



g(z) := g{z) - log^^ - log^ — ^ 

z + c z — b 



2Az^ - log 



z + b 



log 



z — c 



(38) 



' z + c z — b 

and is a contour encircling the three distinct intervals of L. Deforming cz, to a contour 
at infinity and three contours encircling the point z and the two intervals [—b, — c] and [c, 6], 
respectively, one can determine Hs to be 



H,{z) = ~g{z)-2^{z^-a^){z'-b^){z'-c') (a + j' 



\d\ 



A2)^(a2-A2)(62-A2)(A2-c2) 

(39)' 
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and 



XdX 



(^2 _ A2)^(a2 - A2)(62 _ A2)(A2 - c^) 

(40) ' 

Using the definition of Hs{z) in (36), we arrive at the following expansion for Hs{z) at large 
z 

H,{z) = i + 1 r p,(A)A^c^A + ^F'^{A) + • • •, (41) 

Z Z J— a Z 

where Fg is the free energy in the strong region. Therefore, one can expand 
Hs{z)/^J (^2 — a2)(^2 _ _ (o2j g^j^^j demand that it behaves like l/z'^ at large z. This 

gives 

(42) 



XdX 



c ^(a2-A2)(62-A2)(A2-c2) 



and 
A 



-{a' + 6^ + c^) + -(a^fe^ + b^c^ + c\^) 
4 2 



^(a2- A2)(62- A2)(A2-c2) 



1. 



(43) 



Finally, the 1/z^ coefficient of Hs{z) is F^{A). So, 



FVA) = - f -(a« + 6« + c8) - -{a%^ + a^c^ + fe^c^) + {a'b^c^ + a^c^b" + b^c^a^) 
161- 4 2 



+- 



+ 6« + c*^ - (a^6^ + a^c^ + 6^a^ + 6^c^ + c2a^ + c'b'') + 2a%\^ 



b XdX 



+ 4 



where 



+ 6^ + - 2(0^62 + a^c" + b^c") ] ^ 



X^dX 



+ (a' + b' + c' 



^ Jc R(X) Jc 



R{X) 



R{X) := ^(a2- A2)(62_ A2)(A2-c2). 



R(X) 

b X'^dX 

wr 

(44) 
(45) 



Note that (42) and (43) are just two equations for three unknowns a, b and c. This is 
the same situation encountered in [15]. There, the authors have considered an asymmetric 
ansatz for the density of YM2, and have found three equations for four parameters. We will 
discuss this problem in [16]. Now, for the case in hand, the third equation is (35), which 
using eqs. (36) and (37) can be rewritten as 



aJ^ R{z)dz+ I dz P 



b R{z)XdX 
{z^ - A2)i?(A) 



0. 



(46) 
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Now we have three equations ( (42), (43), and (46) ) for three unknowm parameters a, b and 
c. 

To study the structure of the phase transition, we use the following change of variables 

c = s{l-y), 

b^s{l + y), (47) 
a — + e. 

Here, e and y arc equal to zero at the transition point {A = A^). Now, one can expand the 
equations (42), (43) and (46) and obtain 

4^5^ h (As^ H e h 2^5^ H ^ ^ = 0, 48 

s ^ 2s' 8s ^ 4s 8s 64 s ' ^ ' 

lAs' - 1 - TTS + (4As^ + ^)e + (^ - ^)e2 + (lOAs^ - ^)y^ 

2 4 o 4 

+ {As' + ^)ey^ + (2.1s^ - ^7rs)y^ = 0, (49) 

and 

27r - 8^s2 - AAes^ - 27re + ^e^s^ + 2T:e^ + |(16As2 + 177r - 8^es^ + 467re)y2 

491 

+ (2yls^ + — 7r)y^ = 0, (50) 

respectively. These expansions are up to order y^, or (as we will see e is of the order y^). 
Solving eq. (48) for s results 

from which, one obtains ( using eq.(50) ) 

e^-y-"- -y\ (52) 

This shows that e is, in fact, of the order y"^. Substituting these into the eq.(49), we arrive 
at 

.^^f+^.^ (63) 

and 

e = -a H a^, (54) 

6 192 ' ^ ' 

where a is the reduced area, a — ^"j^" ■ Now F'^{A) in eq. (44) can be calculated. The result 

is 

i^i(^) = ^(l + ^«' + ---)- (55) 



Comparing this with FI^{A), eq.(25), it is seen that 

F'siA) - F'^{A) = l^i^^f + (56) 

This shows that we have a third order phase transition, which is the same as the ordinary 
YM2 [10]. 

Acknowledgement M. Ahmohammadi would hke to thank the research vice-chancellor 
of Tehran University, this work was partially supported by them. 



References 

[1] M. Blau and G. Thompson, "Lectures on 2d Gauge Theories", Proceedings of the 1993 
Trieste Summer School on High Energy Physics and Cosmology, World Scientific, Sin- 
gapore (1994) 175-244. 

[2] E. Witten, J. Geom. Phys. 9 (1992) 303. 

[3] M. R. Douglas, K. Li and M. Staudacher, Nucl. Phys. B240 (1994) 140. 

[4] A. Migdal, Zh. Eskp. Teor. Fiz. 69 (1975) 810 (Sov. Phys. JETP 42, 413). 

[5] O. Ganor, J. Sonnenschein, and S. Yankielowicz, Nucl. Phys. B434 (1995) 139. 

[6] M. Khorrami and M. Ahmohammadi, Mod. Phys. Lett. A12 (1997) 2265. 

[7] D. J. Gross, Nucl. Phys. B400 (1993) 161. 

[8] D. J. Gross and W. Taylor, Nucl. Phys. B400 (1993) 181. 

[9] B. Rusakov, Phys. Lett. B303 (1993) 95. 

[10] M. R. Douglas and V. A. Kazakov, Phys. Lett. B319 (1993) 219. 

[11] B. Rusakov and S. Yankielowicz, Phys. Lett. B339 (1994) 258. 

[12] V. S. Popov and A. M. Perelomov, Soviet Math. Dokl. 8 (1967) 712. 

[13] E. Brezin, C. Itzykson, G. Parisi and J. B. Zuber, Commun. Math. Phys. 59 (1978) 35. 

[14] D. Gakhov, "Boundary problems", Russian edition (Nauka, 1975); 

A. C. Pipkin, "A course on integral equations" (Springer, Berhn, 1991). 

9 



[15] J. A. Minahan and A. P. Polychronakos, Nucl. Phys. B422 (1994) 172. 

[16] A. Aghamohammadi, M. Alimohammadi and M. Khorrami, "Uniqueness of the mini- 
mum of the free energy of the 2D Yang-Mills theory at large N" , |hep-th/970708U 



10 



